Carrier dynamics and coherent acoustic phonons in nitride heterostructures 



G. D. Sanders and C. J. Stanton 
Department of Physics, University of Florida, Box 118440 
Gainesville, Florida 32611-8440 
(Dated: February 6, 2008) 

We model generation and propagation of coherent acoustic phonons in piezoelectric InGaN/GaN 
multi-quantum wells embedded in a pin diode structure and compute the time resolved reflectivity 
signal in simulated pump-probe experiments. Carriers are created in the InGaN wells by ultrafast 
pumping below the GaN band gap and the dynamics of the photoexcited carriers is treated in a 
Boltzmann equation framework. Coherent acoustic phonons are generated in the quantum well 
via both deformation potential electron-phonon and piezoelectric electron-phonon interaction with 
photogenerated carriers, with the latter mechanism being the dominant one. Coherent longitudinal 
acoustic phonons propagate into the structure at the sound speed modifying the optical properties 
and giving rise to a giant oscillatory differential reflectivity signal. We demonstrate that coherent 
optical control of the differential reflectivity can be achieved using a delayed control pulse. 

PACS numbers: 63.20.Kr, 63.22.-|-m, 78.20.Hp, 78.47.-|-p 



I. INTRODUCTION 

Femtosecond transient reflectivity spectroscopy has 
proven useful in the study of carrier dynamics in bulk 
semiconductors and semiconductor heterostructures as 
well as the study of the generation and propagation of 
coherent phonons in a number of materials. In partic- 
ular, coherent optical phonons have been observed in 
bulk semiconductors^'^ and coherent acoustic phonons 
have been detected in InGaN/GaN-based semiconductor 
heterostructures,'^— '^•^•^•^ 

We recently developed a detailed theory for carrier 
dynamics and time resolved differential reflectivity in 
t"wo-color pump-probe experiments on nonpiezoelectric 
Ina;Mni_a;As/GaSb heterostructures, -which agreed "well 
■with experimental measurements.^'^" Oscillations "were 
observed in the differential reflectivity which we at- 
tributed to the generation and propagation of coherent 
acoustic phonon wave packets, which altered the local 
dielectric function as they propagated through the struc- 
ture at the LA phonon sound speed. The carrier dy- 
namics in the two-color pump-probe experiments were 
modeled using a Boltzmann equation formalism including 
photogeneration of carriers by the pump laser and their 
subsequent cooling and relaxation by emission of con- 
fined LO phonons. The recombination of electron-hole 
pairs via the Shockley-Read carrier trapping mechanism 
was included in a simple relaxation time approximation. 

The differential reflectivity was obtained by solving 
Maxwell's equations throughout the structure using a 
transfer matrix method. In addition to the coherent 
phonon induced reflectivity oscillations, a strong back- 
ground signal was observed which we attributed to: (1) 
enhanced Drude absorption by photoexcited carriers, (2) 
relaxation dynamics due to cooling of hot photoexcited 
carriers by LO phonon scattering, and (3) nonradiative 
recombination of electron-hole pairs at midgap defects. 

A microscopic theory for the generation and propaga- 
tion of coherent acoustic phonons in piezoelectic wurtzite 



semiconductors was developed in Ref. O (see the erra- 
tum in Ref. and the review article by Chern et al 
m Ref. H. It was found that generation of coherent 
acoustic phonons by means of photoexcited carriers oc- 
curs through the deformation potential electron-phonon 
interaction and the piezoelectric electron-phonon inter- 
action, with the latter typically being an order of magni- 
tude stronger. While the work of Ref. 11 concerned itself 
with the generation and propagation of coherent acous- 
tic phonons in piezoelectric heterostructures, it did not 
address the question of their detection in time resolved 
differential reflectivity experiments. 

Our earlier paper (Ref.0) described the generation and 
propagation of coherent acoustic phonons in Zinc Blende 
semiconductor heterostructures in which the generation 
of coherent acoustic phonons can only occur by means of 
the weaker deformation potential electron-phonon mech- 
anism. In this paper we are motivated to extend the re- 
sults of Ref. ,9, to piezoelectric wurtzite heterostructures 
where the coherent phonon induced differential reflectiv- 
ity oscillations are expected to be much stronger. 



II. THEORY 

We model photogeneration of electrons and holes and 
the subsequent excitation of coherent acoustic phonons 
in a nitride based multi-quantum well (MQW) pin diode 
shown schematically in Fig. ^ Such a nitride structures 
have been used in recent experimental studies by Y. D. 
Jho et al in Refs. ITsI andiM The structures are grown 
on an undoped 1.5 /Ltm thick GaN layer grown on top of 
an flc-plane sapphire substrate and we base our theoret- 
ical studies on a similar structure. The intrinsic region 
of the diode is an undoped MQW heterostructure con- 
sisting of five 22 A Ino.15Gao.85N wells separated by 100 
A GaN barriers. The growth direction is along (0001). 
The multi-quantum well structure is the intrinsic region 
of a pin diode with contacts consisting of p- and n-doped 
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FIG. 1: Schematic diagram of the multi-quantum well pin 
diode structure. 



GaN substrate and cap layers. The thickness of the n- 
GaN cap layer is taken to be 100 A, which is thinner 
than that used in Refs. andll^. and the width of the 
intrinsic region L is 710 A as can be seen in Fig. ^ The 
voltage drop across the intrinsic re gion is AT^. For the 
structure described in Refs. and ll4l the voltage drop 
across the intrinsic region is At^ w 1.6 volts when the 
external voltage across the diode is 1^ = volts and we 
adopt this value in our studies. 



The Hamiltonian for the valence bands is a 6x 6 matrix. 
Following Ref. 17, the Hamiltonian (relative to the top 
of the valence band) can be block diagonalized into two 
degenerate 3x3 submatrices if we adopt the Bloch basis 
states 

\v,l) = -^\{X + ^Y)]) + ^\{X-^Y)i) (3a) 



^\[X-^Y)])^^\{X + ^Y)i) (3b) 



\v,i) = p*\Z])+p\Zi) 
\vA) 



(3c) 



l«,6) 



-^\{X + ^Y) t)-^l(X-zy) i) (3d) 
^\{X-^Y)]) + ^\{X + ^Y)i) (3e) 



= ~P*\Z^)+(i\Z i). 



(3f) 



The phase factors, a and (3, are functions of the angle 
(j) = tan^^ {ky / kx) and are given by 
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The block diagonalized Hamiltonian can be written as 



Bulk Bandstructure 



where the upper and lower blocks of the Hamiltonian are 



Our first task is to model the electronic states in the 
device shown in Fig. ^ In bulk systems, the conduction 
and valence bands in wurzite crystals including the effects 
of strain are treated using effective mass theory. Near the 
band edge, the effective mass Hamiltonian for electrons 
is described by a 2 x 2 matrix which depends explicitly 
on electron wavevector, k, and the strain tensor, e. The 
electron Bloch basis states are 
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and 



F Kt -iHt 
Kt G ^~iHt\ (6) 
iHt A + iHt A 



F Kt iHt 
H^^,{k,e)=\ Kt G A + iHt 
-iHt A -iHt A 



(7) 



(1^) The elements appearing in the 3x3 Hamiltonian ma- 

(Ib) trices are 



The conduction band effective mass Hamiltonian is di- 
agonal and is given by (relative to the bottom of the 
conduction band) fRefs. Iisl and lT^ 



^f2x2(k,e) = { 



2j,2 



h-'k 



2m: 



2m* 



xy 

yy)} I2x2- 



(2) 



where I2x2 is the identity matrix. The electron effective 
masses along z (taken parallel to the c-axis) and in the 
xy plane are m* and m*j^, respectively, fc^ = 



and a, 



yy and e^z 



and 

are strain tensor components, and ac.z 
c.xy are the deformation potentials. 
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A = \/2 Aa 
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TABLE I: Material parameters for wurtzite GaN and InN. 
Values are taken from Ref. 18 unless indicated otherwise. 
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Here, mo is the free electron mass, the A'^s are efFective 
mass parameters, the D[s are the Bir-Pikus deformation 
potentials, and the A's are related to the crystal field 
splitting, Acr, and spin-orbit splitting, Aso, by Ai = Acr 
and A2 = A3 = A,o/3. 



B. Material parameters 

The material parameters for GaN and InN are shown 
in Table m Unless otherwise stated, the parameters are 
taken from Ref. 

The temperature dependent direct band gaps are given 
by the empirical Varshni formula^ 



EgiT) 



(15) 



where Eg is the band gap at absolute zero, T is the 
temperature and ag and (3g are the Varshni parame- 
ters. These band gap parameters are listed in Table ^ 
for wurzite GaN and InN. 

For the In^jGai-^^N alloy the temperature and compo- 
sition dependent band gap is given by 

Eg(x,T) = xEgjnN + (1 - x)Eg^GaN - bx{l ~ x) (16) 

where EgjnN and Eg^QaN are the temperature depen- 
dent band gaps in bulk InN and GaN as given in Eq. [T51 
and 6 = 3.0 eV is the band gap bowing parameter for 
Ina;Gai_a;N.^* While the value of the GaN band gap is 
well known, there is some uncertainty in the value of 
the InN band gap used in computing the band gap in 
Ina;Gai_a;N. The InN band gap determined from absorp- 
tion measurements on InN samples range from 1.7 to 3.1 
eV with most values centered around 1.9 eVr^ Based on 
optical measurements on ultrapure InN samples, it has 
recently been suggested2Si2Si2L2& that InN is a narrow 
gap semiconductor with a band gap of 0.7 eV as opposed 
to 1.994 eV as cited in Ref.IlM Nevertheless, for the small 
values of the Indium concentration, x, that we use in this 
study, the computed band gap for Ina;Gai_2;N should be 
rather insensitive to the assumed value of the InN band 
gap and we adopt the InN band gap and band gap bowing 
parameter for the Ina;Gai_a;N alloy cited in Ref.[l^ Since 
the band bowing parameter, 6, for Inj;Gai_a:N based on 
an InN band gap of 0.7 eV has not been redetermined. 



Parameter 



Mass Density 
Po (gm/cm^) 
Lattice constants 
ao( A) 
co( A) 

Direct band gaps 
E, (eV) 
ag (meV/K) 
(K) 

Hole splitting energies 
A„. (eV) 
(eV) 

Electron effective masses 
m* (mo) 
mly (mo) 

Hole effective mass parameters 

A2 
Ai 

A„ 
Ae 

Electron deformation potentials 
ac,z (eV) 
ac,xy (eV) 

Hole deformation potentials 
Di (eV) 
D2 (eV) 
Di (eV) 
Di (eV) 

Elastic stiffness constants 
Cii (GPa) 
C12 (GPa) 
Ci3 (GPa) 
C33 (GPa) 

Piezoelectric constants 
ei3 (C/m^) 
633 (C/m^) 
Dielectric constants 
£0 
£ 00 

Refractive index 
n 



"Ref. 



''Ref. 2i 

=Ref. 21 

''Ref. 25 

■^Ref. 23 



1£ 



GaN 



6.15 

3.189 
5.185 

3.507 
0.909 
830 

0.019 
0.014 

0.2 
0.2 

-6.56 
-0.91 
5.65 
-2.83 
-3.13 
-4.86 

-6.5 
-11.8 

-3.0 
3.6 

8.82 
-4.41 

390 
145 
106 
398 

-0.35 
1.27 

8.9 " 
5.7 " 

2.27 



InN 



6.81 

3.545 
5.703 

1.994 
0.245 
624 

0.041 
0.001 

0.12 
0.12 

-8.21 
-0.68 
7.57 
-5.23 
-5.11 
-5.96 

-3.5 
-3.5 

-3.0 
3.6 
8.82 
-4.41 

223 
115 
92 
224 

-0.57 
0.97 

15.3 ' 
8.4 " 

2.71 ' 
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our use of the older parameters has the virtue of including 
band bowing effects in our calculations for small values 
of X. 

We obtain electron effective masses for IwxG&i-x^ by 
linearly interpolating the reciprocals of the masses as a 
function of x, i.e., the concentration dependent effective 
masses are taken to be 
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For all other material parameters in TableQ] we use linear 
interpolation in the composition to obtain values for the 
alloy. 



C. Carrier states in multi-quantum well diode 

In computing the electronic states in the intrinsic re- 
gion of the pin diode shown in Fig.^ we take the barriers 
in the ohmic contacts to be infinite for simplicity. The 
quantum confinement potentials for electrons and holes 
in the intrinsic region arise from bandgap discontinuities 
between well and barrier regions and the strain-induced 
piezoelectric field. Thus the confinement potential is 



Va{z,t) = VQ,gap(2:) -I- Vp\c7.o{z), 



(18) 



where a = {c, t;} distinguishes between electrons and 
holes. 

If the position-dependent band gap in the MQW is 
Eg{z) as determined from Eg. 1151 then the confinement 
potentials for electrons and holes are given by 



^,gap(^) 



E. 



5, mm 



{Eg{z) - Eg^,-^^) (19a) 

K,gap(2) = -(l-Qc) {Eg{z) ~ Eg^^,,^) (19b) 

where -Eg^min — niinz[i?g(z)] is the minimum of the posi- 
tion dependent band gap and — 0.6 is the conduction 
band offset iSS 

The confinement potential due to the strain-induced 
piezoelectric field satisfies 



''^P'"^"^^^ = \e\ El{z), 



(20) 



where |e| is the electric charge and El{z) is the strain- 
induced piezoelectric field. The piezoelectric field in the 
diode is obtained from the requirement that the electric 
displacement vanishes,— Thus 



E'Az) 



47r 

£0(2) 



(21) 



where Pq is a constant polarization induced by the ex- 
ternally applied voltage, AT^, and eo(z) is the position- 
dependent static dielectric constant. The value of Pq 



is obtained from the voltage drop across the diode (of 
length L). The voltage drop in Fig. ^ between source 
and drain due to the induced piezoelectric field is just 



Ay = - 



dz E°{z), 



(22) 



from which Pq can be determined. The magnitude of the 
strain-induced polarization directed along the z-direction 
is given by 

P°(z) = 631 (z) (e^^(z) + eyy{z)) + 633(2)6^^(2), (23) 

where 631(2) and 633(2) are the position-dependent piezo- 
electric constants and exxiz), ^yy(z) and tzz(z) are the 
position-dependent strain tensor components. The orien- 
tation of the strain- induced polarization, P^[z), is such 
that the piezoelectric field opposes the built-in electric 
field in the intrinsic region as has been discussed in sev- 
eral previous studiesi^^i^^i^^i'^^i^'^i'^^ 

In a pseudomorphically strained MQW diode, the GaN 
source and drain contacts are assumed to be unstrained 
while the in-plane lattice constants in the In^^Gai-^jN 
wells and GaN barriers adjust to the lattice constant in 
the source and drain. For a MQW grown along the [0001] 
direction (the z-direction) the z-dependent strain is given 



^Xx{z) — €yy{z) 



flo — a(z) 
a[z) 



(24) 



Here ao is the lattice constant in the GaN source and 
drain and a{z) is the z-dependent lattice constant in the 
MQW active region. Minimizing the overall strain, we 
find^a 



2 Ci3(2) 
C33{z) 



(25) 



where 6*13(2) and 6*33(2) are z-dependent elastic stiffness 
constants. 

The diode structure breaks translational symmetry 
along the 2 direction. Thus, quantum confinement of 
carriers in the MQW active region gives rise to a set of 
two-dimensional subbands. The effective mass wavefunc- 
tions are 



A 



(26) 



where a = {c, v} refers to conduction or valence sub- 
bands, n is the subband index, k — {kx, ky,Q) = (fc, (j>) is 
the two-dimensional wavevector, and j labels the spinor 
component. For conduction subbands, (a=c) j = 1,2 
while for valence subbands (a=v) j = 1...6. The slowly 
varying envelope functions F^j^ -{z) are real and depend 
only on fc = |k|, while the rapidly- varying Bloch basis 
states |a, j) are defined in Eqs. Q and The valence 
Bloch basis states depend on the orientation, (p, of k in 
the xy plane as given in Eq. Q. The area of the MQW 
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sample in the xy plane is A, and p = {x, y, 0) is the pro- 
jection of r in the plane. 

The envelope functions satisfy a set of effective-mass 
Schrodinger equations 

J2 {HfAk) -f 5,^,, [14 (z) - E^{k)]} Fl^^^,{z) = 0, 

3,0' 

(27) 

subject to the boundary conditions 



(28) 



where Va{z) are the quantum confinement potentials for 
conduction and valence electrons defined in Eq. 1181 and 
Ef^{k) are the energy eigenvalues for the nth conduction 
or valence subband. Note that in the envelope function 
approximation, the subband energy depends only on the 
magnitude k of the transverse wavevector and not on the 
angle cf). The matrix operators H°'-, (k) depend on z and 

are obtained by making the replacement k^ a-nd 
letting all material parameters be z-dependent operators 
in the matrices H''{k,e) given in Eqs. and To 
ensure the Hermitian property of the Hamiltonian, we 
make the operator replacements3i 



d'^ d d 



and 



dz 



dz' 



oz oz 



(29) 



(30) 



We arrive at a set of coupled ordinary differential equa- 
tions (ODE's) subject to the two-point boundary value 
condition of Eq. H28|) . These are solved for the envelope 
functions and subband energies. In practice, we intro- 
duce a uniform grid, {zi}, along the z-direction and finite- 
difference the effective mass Schrodinger equations to ob- 
tain a matrix eigenvalue problem which can be solved 
using standard matrix eigenvalue routines. The result- 
ing eigenvalues are the subband energies, E^{k), and the 
corresponding eigenvectors are the envelope functions, 
F"i^j{zi), defined on the finite difference mesh. 



electron in subband n with wave vector k. The Boltz- 
mann equation including photoexcitation of hot electron- 
hole pairs by the pump, the subsequent cooling of these 
carriers by emission of confined LO phonons, and the 
recombination of electron-hole pairs is 



-/,?(k)Kk' [l-r'(k')]} 



a/,?(k) 



dt 



(31) 



To simplify the calculations, we use an axial approxi- 
mation in which the distribution functions are replaced 
by their angular averages in the x-y plane of the multi- 
quantum well^*^ Thus the time dependent conduction 
and valence electron distribution functions, f^{k) and 
f^,{k), are functions of the magnitude of k. 

The first two terms on the right hand side of Eq. I|31|) 
describes the rapid cooling of hot electrons and holes by 
emission and absorption of confined LO phonons. Rapid 
LO phonon scattering is an intraband process which cre- 
ates quasi-thermal electron and hole distributions at their 
respective band edges. The confined LO phonon scat- 
tering rate in the multi-quantum well, W^'^, , given in 
Refs.lHandi, is the rate at which a conduction (or va- 
lence) electron in subband n with wave vector k scatters 
to subband n' with wave vector k'. 

The last term on the right hand side of Eq. H31|l consists 
of two terms 



p/„"(k)l 






-1- 


r9/„"(k)i 








pump 





(32) 



rela 



and describes the change in the conduction or valence 
electron distribution function due to the action of the 
pump as well as carrier relaxation through recombination 
of electron-hole pairs. We assume that electron- hole pairs 
can recombine with a phenomenological time constant tq. 
Expressions for the photogeneration rate can be found in 
Ref. 0- It is assumed that the pump laser is a Gaussian 
pulse with an intensity FWHM of Tp peaked at t = with 
a Lorentzian spectral lineshape with a FWHM of 7. 



D. Carrier dynamics 



E. Optical properties 



In two-color time resolved differential reflectivity ex- 
periments a pump laser is used to excite electrons from 
the valence to the conduction subbands of the multi- 
quantum well. The photoexcited carriers then relax 
through scattering, changing the optical properties of 
the heterostructure in the process. We simulate these 
processes using Boltzmann transport equations which we 
solve numerically For each electron and hole subband 
state with energy, E!j^{k), we have a time dependent dis- 
tribution function, /"(k, i), which gives the probability, 
as a function of time, of finding a conduction or valence 



For dipole transitions between conduction and valence 
electronic states near the band edge, we calculate optical 
properties using Fermi's golden rule. The contributions 
to the real and imaginary parts of the dielectric function 
from the effective mass band edge states are given by 



ei{fujj) = 



8ne- 



^ J2 |e.<x(k) ir^,ik)-rjk)) 

n,n' ,k 



(33) 



6 



and 



V 



n,n' ,k 



X <5(A£;:-;,(fc)-M 



(34) 



where AE^^,^,{k) = i;fj(fc) - i;;j,(fc) are the fc-dependent 
transition energies between conduction and valence sub- 
band states, V = A L is the sample volume, and e is the 
unit complex polarization vector. Expressions for the 
dipole matrix elements d^'"^, (k) for dipole allowed tran- 
sitions between conduction subband n and and valence 
subband n' can be found in Ref. |0- 

In Eqs. and the spectral lineshapes are de- 
scribed by Dirac delta functions. In practice we replace 
the delta functions with Lorentzian lineshapes having a 
full width at half maximum (FWHM) of F. The real and 
imaginary parts of the dielectric function, including line- 
shape broadening, can be obtained from Eqs. (|33|l and 
using the operator replacements'^'' 



and 



1 



r/(27r) 



AE - hw 



AE-huj {AE - hujy + (T/2y 



(35) 



(36) 



F. Generation and propagation of coherent 
acoustic phonons 

The ultrafast photogeneration of electrons and holes 
in the InGaN multi-quantum wells by the pump gives 
rise to coherent longitudinal acoustic (LA) phonons 
which propagate into the diode structure. Coherent 
acoustic phonons give rise to a macroscopic lattice 
displacement i^i^i^iSiLSiiiiiS Since the photogenerated car- 
rier distributions are functions of z, the transient lattice 
displacement U{z, t) due to photogenerated carriers is in- 
dependent of X and y and is parallel to z. As discussed 
in Refs. 13 and|3, U{z,t) satisfies a loaded string equa- 
tion. In the presence of a position dependent longitudinal 
acoustic sound velocity, Cs{z), we have 



d^U{z,t) d 



dU(z,t) 
dz 



Siz,t) (37) 



where S{z, t) is a driving or loading function and depends 
on the photogenerated carrier distributions. The longi- 
tudinal acoustic sound velocity is given by 



Cs{z) 



Po{z) 



(38) 



where 6*33(2) and Pq{z) are the position dependent elastic 
stiffness constant and mass density. 



The loaded string equation is solved subject to the ini- 
tial conditions 



U{z,t = -co) 



dU{z,t 



'OO 



dt 



0. 



(39) 



We solve the loaded string equation numerically by finite 
differencing Eq. H37|l inside a computational box whose 
left edge, z^, is the semiconductor-air interface in the 
p-GaN contact and whose right edge, zr, lies inside the 
n-GaN substrate (see Fig. ^ . At z/j we impose absorbing 
boundary conditions while at there are no perpendic- 
ular forces at the semiconductor-air interface. Thus we 
solve the initial value problem subject to the left and 
right boundary conditions 



dU{zL,t) 



dz 



= 



and 



dU{zR,t) 
dz 



Cs{zr) 



dU{zR,t) 
dt 



0. 



(40a) 



(40b) 



Starting with the second quantized Hamiltonian for the 
electron-phonon interaction, a microscopic expression for 
the driving function was derived in Ref. JTi using the den- 
sity matrix formalism (see also the erratum in Ref. as 
well as the review article in Ref. In wurzite materials 
such as InGaN the electron-phonon interaction contains 
contributions from piezoelectric and deformation poten- 
tial coupling. Under typical experimental conditions the 
microscopic expression for the driving function can be 
simplified toii*i^ 



S{z,t)=Y,S,{z,t) 



(41) 



where the summation index, z/, runs over carrier species, 
electrons, heavy holes, light holes and split off holes. 



i.e 



The partial driving functions including both piezoelectric 
and deformation potential coupling are given by 



I f d 

Si,{z,t) = ±— a^ — 
Po V dz 



47r|e| 633 



p,{z,t) (42) 



where po is the mass density, 633 is the piezoelectric con- 
stant, is the high frequency dielectric constant and 
Pi,{z,t) is the photogenerated carrier density. The de- 
formation potentials are a^. For conduction electrons, 
— 0'c,z, for heavy and light holes, = Di + D3, 
and for crystal field split holes, a,^ = Di. We note that 
Eq. 142|) was derived in the elastic continuum limit by 
Chern et al^ The driving function satisfies the sum rule 



dz S{z,t) = 0. 



(43) 



as shown in Refs. OandU 

Since the photoexcited holes are predominantly a mix- 
ture of heavy and light holes, we use a^, = D13 = D1+D3 
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in computing hole deformation potential contributions in 
Eq. H41|l . The sum over species, then yields the total 
driving function 



S(z,t) = — <ac,z^ 



47r|e|e33 



1 

Po 



D 



13 



d_ _^ 47r|e|e33 
dz e„o 



Pciz,t) 



Ph{z,t), (44) 



where pc{z,t) and p]^[z^t) are the photogenerated con- 
duction electron and valence hole densities. 

The photogenerated electron and hole densities are ob- 
tained from the envelope functions defined in Eq. H26(l 
and the time dependent electron and hole distribution 
functions, /^(fc,t) and /j",(fc,t), obtained from solving 
the coupled Boltzmann equations We have 



Pe{z,t) 



(45a) 



A 



and 



Phiz, t)^ J E (/„"?(fc)-/:'(fc,i))|i^,"',fe,,W|'- 

7i' ,j,k 

(45b) 

where f^ik) and f^{k) are the initial Fermi-Dirac dis- 
tribution functions for conduction and valence band elec- 
trons in subbands n and n' . 

The propagating coherent phonon displacement field 
U{z,t) gives rise to a propagating strain with 



Az,t) 



dU{z,t) 
dz 



(46) 



The propagating strain field alters the optical properties 
of the sample which is then detected by the delayed probe 
pulse. 



G. Transient probe response 

To compute the time dependent probe transmission 
and reflection coefficients we need to model the dielec- 
tric function in the diode structure shown in Fig. ^ We 
will denote this time and position dependent dielectric 
function as 

e{huj, z, t) = ei{hLL>, z,t) + i £2(fujj, z, t) (47) 

where hiu is the probe energy. We can solve Maxwell's 
equations for the time dependent probe reflection coeffi- 
cient, R^hio, t), usin g th e transfer matrix method of Ref.'^ 
as described in Ref. |33. 

There are several processes which contribute to the 
dielectric function in the InGaN/GaN multi-quantum 
wells. The first contribution to the dielectric function 
is a Drude term due to photoexcited free carriers which 
gives a contribution to the dielectric function of 



e{huj, z, t)i 



{hujp{t)Y 
{hu;f 



(48) 




3 4 5 6 7 
Photon Energy (eV) 



FIG. 2: (Color Online) Model dielectric function for bulk GaN 
at T = K as a function of photon energy used in calculating 
the dielectric function in the diode source and drain contact 
regions. 



where ujp is the plasma frequency. The Drude contri- 
bution to the dielectric function in Eq. H48|l is uniform 
in the multi-quantum well structure and vanishes every- 
where else. In the random phase approximation (RPA), 
the time dependent plasma frequency is given bjiS^ 



47re^ 



n,k 



m, 



E 

n' ,k 



ml,{k) 



(49) 



where L = 710 A is the width of the multi-quantum well 
structure in Fig. ^ and A is the cross sectional area of 
the diode. The effective masses for electrons and holes 
are obtained from the computed energy bands as 



1 



(50) 



Note, that in our convention, the bandstructure is com- 
puted in the electron picture so that the hole masses 
in Eq. H50|l are usually negative. Thus both terms in 
Eq. (|49|l give positive contributions to the squared plasma 
frequency. 

Another contribution to the multi-quantum well di- 
electric function, which we denote e{hu!, z, t) mqw, comes 
from interband transition between the effective mass 
band edge states. The real and imaginary parts of 
e{hio, z,t)MQw are given in Eqs. (|33|) and (|34|l and are 
time dependent due to variations in the electron and hole 
distribution functions with time. 

There is also a background dielectric function, Sbihuj) 
due to all the higher lying electronic transitions. For 
simplicity, we treat these contributions to the dielec- 
tric function in the multi-quantum well region using the 
model dielectric function for bulk GaN in Ref. |2 with 
contributions from the Eq and Eq + Aq critical points 
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removed. These correspond to contributions from the 
muhi-quantum weU effective mass electronic states and 
are aheady included in Eqs. and 

We assume the coherent phonon alter the background 
dielectric function because of strain induced variations in 
the energy gaps associated with each transition^i*^ Den- 
sity functional calculations of the deformation potentials 
for the El transitions in a number of semiconductors^ 
have shown that the deformation potentials associated 
with the Eq and Ei features are equal to within 20%. 
Thus, at photon energies at or below the Eq and Ei 
transitions we assume that the effect of temperature and 
strain on Sb is to introduce a rigid shift in the dielectric 
function such that 

ebihuj,z,t) = Ebihw - AEg{T) - acviexx + £yy + £zz))- 

Here, is an effective deformation potential which we 
take to be acv — —18.3 eV. The coherent phonon strain 
tensor components are, e^x, £yy, and e^z, and AEg(T) = 
Eg{T) — Eg is the band gap shift due to temperature 
variations with Eg{T) being the temperature dependent 
band gap defined by the empirical Varshni formula in 
Eq. ini). 

The total dielectric function in the multi-quantum well 
region is obtained by adding the Drude, multi-quantum 
well, and background contributions in Eas. (jI5)l . (O, 
(Oljl and (EH), i.e. 

e{hu, z, t) — e{hiu, z,t)D + e{huj, z, t) mqw + z, t) 

(52) 

In the GaN contact regions, we use the model dielectric 
functions for bulk GaN described in Ref. Fig. |21 

shows the real and imaginary parts of the model dielectric 
function for bulk GaN at T = K in the absence of 
strain. Temperature and strain effects are included using 
the same rigid shift model as defined in Eq. (|51|l . Note 
that the dielectric functions in the GaN source and drain 
regions are modulated by the coherent phonon strain field 
as it propagates through the structure. 

III. RESULTS 

A. Multi-quantum well band structure and 
electronic states 

The band diagram for the conduction and valence sub- 
band states in the intrinsic region of our pseudomorphi- 
cally strained InGaN/GaN MQW pin diode structure, 
shown schematically in Fig. can be seen in Fig.|2|as a 
function of position along the growth direction, z. The 
voltage drop across the intrinsic region of the diode is as- 
sumed to be Ay w 1.6 Volts which is the experimentally 
measured value obtained in similar structures described 
in Refs. and ^3 for an externally applied voltage of 
V Q Volts. The center of the MQW intrinsic region 
is defined to be at z = and the valence band poten- 
tial is set to zero at this point. Since the in-plane lat- 




-400 -200 200 400 

Position (Angstroms) 

FIG. 3: (Color Online) Band diagram for conduction and 
valence bands in the intrinsic region of an InGaN/GaN MQW 
pin diode structure. For an external bias of 1/ ~ Volts, the 
voltage drop across the intrinsic region of the diode is assumed 
to be Al/ ~ 1.6 Volts as described in the text. 



tice constant throughout the structure is the same as in 
the GaN source and drain contacts, the GaN barriers in 
the MQW region are unstrained while the Ino.15Gao.85N 
wells are under biaxial compression. The electric field is 
2.118 Mv/cm in the InGaN wells and -0.686 Mv/cm in 
the GaN barriers. 

Using the confinement potentials in Fig. |31 we can 
solve the effective mass Schrodinger equation (|77|) for the 
electronic subbands E"{k) and the envelope functions, 
^nkji^)- We can then compute the optical dipole ma- 
trix elements d^'"^, (k) . The conduction and valence sub- 
bands, E"{k), for the confinement potentials in Fig.|3|are 
shown in Fig. ^ as functions of k in the axial approxima- 
tion. Note that the scale is different for the conduction 
and valance band levels in the upper and lower panels. 



B. Multi-quantum well differential absorption 

Given the conduction and valence band distribution 
functions /^(fc) and /^(fc), the real and imaginary parts 
of the dielectric function as a function of photon energy 
can be obtained from Eas. 1331 and 1341 The room tem- 
perature absorption coefficient in the pin diode MQW 
intrinsic region due to intersubband optical transitions 
between the effective mass states is shown in Fig. |S1 

To illustrate the effects of carrier dynamics on the op- 
tical properties, we simulate a two-color pump-probe ex- 
periment on the diode structure with a pump laser hav- 
ing a photon energy of 3.2 eV (A = 387.5 nm) as indi- 
cated in Fig. 13 The pump fluence is assumed to be .001 
J /cm^ and the FWHM temporal and spectral widths are 
taken to be Tp = 1 ps and 7 = 4 meV respectively. The 
band gap for GaN at 300 K is 3.43 eV and the band gap 
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k (nm"') 

FIG. 4: (Color Online) Conduction and valence subbands for 
the InCaN/CaN MQW pin diode structure of Fig. Q] with an 
external bias oi V ~ Volts. 



for Ino.15Gao.85N is 2.87 eV. We chose the photon en- 
ergy of the pump laser to lie above the band gap of the 
Ino.15Gao.85N wells but below that of the GaN barriers in 
order to preferentially photoexcite carriers in the wells. 

The pump laser creates electron hole pairs in the quan- 
tum wells which then rapidly cool through confined op- 
tical phonon scattering and slowly recombine across the 
band gap. In our simulations we solve the set of Boltz- 
mann equations in Eq. H^lfl for the time-dependent con- 
duction and valence electron distribution functions /^(fc) 
and fn{k). The phenomenological electron- hole recom- 
bination time is take to be tq = 100 ps since, in an 
earlier experimental studyjSi electrons and holes in an 
Ina;Mni_2.As/GaSb heterostructure were found to recom- 
bine at midgap defects on a time scale of tq w 200 ps. The 
time-dependence of the dielectric functions ei{huj) and 
£2(^.0;) is determined by the time-dependence of the con- 
duction and valence band distribution functions, /^(fc) 
and /^(fc). 

The differential absorption coefficient, defined as the 
change in the absorption coefficient with respect to its 
initial thermal equilibrium value, is shown in Fig. as a 
function of delay time and probe photon energy. There 
is a sharp dip in the differential absorption at the pump 
photon energy (3.2 eV) which can be attributed to bleach- 
ing of the direct transition between photogenerated pairs 
of electrons and holes. In the quantum well structure 
with strong built-in piezoelectric fields, a given quantum 
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FIG. 5: (Color Online) Initial absorption coefficient at room 
temperature due to optical transitions between band edge ef- 
fective mass states in an InGaN/GaN MQW pin diode struc- 
ture with an external bias of V" ^ Volts. 




FIG. 6: (Color Online) Time dependent differential absorp- 
tion as a function of delay time and probe photon energy in 
a simulated pump- probe experiment in a InGaN/GaN MQW 
pin diode structure with an external bias of V Volts. 



well state is optically coupled to quantum well states in 
adjacent wells and in the GaN continuum. When the 
pump creates electron-hole pairs, these interband transi- 
tions are also blocked giving rise to the negative differen- 
tial absorption seen at higher photon energies in Fig. 
The relaxation of the differential absorption as a func- 
tion of delay time occurs on a time scale of hundreds of 
picoseconds and is controlled by intersubband recombi- 
nation of the photogenerated electrons and holes. 
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C. Coherent phonon generation and propagation 

In our simulated pump-probe experiment, electrons 
and holes are photoexcited in the conduction and va- 
lence bands by the pump laser. These carriers in turn 
generate coherent acoustic phonons. From the photoex- 
cited electron and hole distribution functions Pe{z^ t) and 
Ph{z,t) defined in Eq. (|45|1 . we obtain the driving func- 
tion <S'(z, t) from Eq. H44|) and solve the coherent phonon 
loaded string equation numerically on a finite difference 
space-time grid. 

The coherent phonon lattice displacement, U{z,t), ob- 
tained from the loaded string equation H37|) gives rise to 
a propagating strain field. The coherent phonon strain 
tensor component ezz{z,t) is shown in Fig. |7|as a func- 
tion of position in the sample for several equally spaced 
delay times ranging from to 100 ps. 

Following photogeneration of carriers by the pump, a 
localized strain appears in the multi-quantum well region 
as can be seen in Fig. [7| This is due to near steady-state 
loading by the driving function at long times. Assuming 
the driving function, S{z,t), is approximately time inde- 
pendent at long times, the loaded string equation l(H7|) 
can be integrated once in the steady-state limit. The 
resulting steady-state strain is 



dUjz) 
dz 



dz' 



S{z') 



(53) 



where Cs{z') is the longitudinal acoustic sound speed 
in the InGaN quantum wells and S{z') is the approxi- 
mately time independent driving function left behind in 
the multi-quantum well at long times. The fact that the 
steady-state strain is localized in the InGaN/GaN multi- 
quantum wells follows directly from the sum rule (|43|l . 

In addition to the localized strain in the InGaN/GaN 
multi-quantum well region of the diode, transient strain 
pulses are seen to radiate into the GaN layers at the 
longitudinal acoustic sound speed. Two transient strain 
pulses are generated, one propagating to the left and the 
other to the right. The leftward propagating pulse is to- 
tally reflected off the semiconductor-air interface in the 
p-GaN layer and trails the rightward propagating pulse 
as it propagates into the n-GaN layer and the GaN sub- 
strate. 



D. Differential reflectivity 

The computed time resolved differential reflectivity in 
our simulation as a function of probe delay is shown in 
Fig. |Slfor probe energies ranging from 3.0 to 4.0 eV in 
increments of 0.1 eV. 

The oscillation observed in the differential reflectiv- 
ity can be attributed to propagation of the strain pulse 
through the diode. The propagating strain tensor shown 
in Fig. [3 alters the local dielectric function as it propa- 
gates. The change in the complex dielectric function is 
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FIG. 7: (Color Online) The strain tensor component as a 
function of position and delay time. 
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FIG. 8: (Color Online) Computed time resolved differential 
reflectivity as a function of probe delay for simulated pump- 
probe experiments on the InGaN/GaN pin diode shown in 
Fig.Q The pump photon energy is 3.2 eV and differential re- 
flectivity is computed for probe-photon energies ranging from 
3.0 to 4.0 eV. The curves are offset for clarity. 
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FIG. 9: (Color Online) Derivative of the complex GaN di- 
electric function with respect to strain as a function of thne 
probe photon energy. The solid line is the real part and the 
dashed line is the imaginary part. 



given by 



Ae{hLLj, z, t) 



de{huj) 



(54) 



The total derivative of the complex dielectric function 
with respect to strain in Eq. H54|l measures the sensitivity 
of the local dielectric function to small changes in the 
strain induced by the coherent phonons and is obtained 
from the GaN model dielectric function by differentiating 
with respect to e^z taking care to eliminate e^x and Syy in 
favor of Szz- The real and imaginary parts of ds/dszz in 
GaN are plotted as a function of the probe photon energy 
in Fig. 1^1 The sensitivity of the dielectric function to 
changes in the local strain depends strongly on the probe 
photon energy and is especially large near the GaN band 
gap. 

The differential reflectivity curves in Fig. |H1 are ob- 
tained by solving Maxwell's equations in the structure 
for the dielectric function in Eq. H47|l . We note that the 
differential reflectivity signal in Fig. [S] is dominated by 
strong coherent phonon oscillations. In an earlier exper- 
imental and theoretical study of time resolved reflectiv- 
ity in Ina;Mni_a;As/GaSb heterostructures (Ref. we 
found similar oscillations in the differential reflectivity 
signal which were a weak perturbation on the background 
signal. The background signal in the earlier study was 
attributed to: (1) enhanced Drude absorption resulting 
from the increase in carriers and (2) the relaxation dy- 
namics associated with the relaxation of highly nonequi- 
librium photoexcited carriers. In the time resolved differ- 
ential reflectivity studies of two-color pump-probe exper- 
iments in InMnAs/GaSb heterostructures described in 
Ref. 1^, coherent acoustic phonons were generated exclu- 
sively through the deformation potential electron phonon 
interaction.^ The strong differential reflectivity oscilla- 



tions seen in the present study, in contrast to the earlier 
one, can be attributed to an order of magnitude enhance- 
ment in the coherent phonon driving function, S{z,t), 
due to the strong piezoelectric electron-phonon interac- 
tion term in Eq. 

The reflectivity oscillations can be qualitatively un- 
derstood as follows. The propagating strain pulse in 
Fig. [3 gives rise to a perturbation in the GaN dielectric 
function which propagates at the acoustic sound speed. 
Our diode structure with the propagating perturbation 
in the dielectric function acts as a Fabry-Perot interfer- 
ometer and the period of the reflectivity oscillations due 
to the propagating coherent acoustic phonon wavepacket 
is approximately^ 



T 



A 



2 Cs n{\) 



(55) 



where A = 27r c/w is the probe wavelength, hu; is the 
probe photon energy, Cg is the LA sound speed in GaN 
and n{X) is the wavelength dependent refractive index. 
The refractive index can be obtained from the GaN model 
dielectric function in Fig.|21as^'* 



"(A) - (ei(A) + y^(A)2+£2(A)^ 



(56) 



The index of refraction in GaN obtained from the 
model dielectric function in Fig. [21 and using Eq. H5t)|) is 
shown in the upper panel of Fig. llOl as a function of probe 
photon energy and the period of the coherent phonon in- 
duced reflectivity oscillations based on the estimate of 
Eq. H55|l is shown as a solid line in the lower panel. The 
period of the reflectivity oscillations obtained from our 
simulations and presented in Fig. |H1 are plotted in the 
lower panel of Fig. ^1 as dots. As we can see, the simu- 
lation results agree very closely with the simple estimate 
in Eq. 

In Fig. |S1 the coherent phonon induced reflectivity os- 
cillations are observed to decay as a function of probe 
delay and we note that the temporal decay of the reflec- 
tivity signal is more rapid at higher probe photon ener- 
gies. This has nothing to do with decay of the propagat- 
ing coherent phonon strain pulse, but rather can be at- 
tributed to the fact that the absorption coefficient in GaN 
is rapidly increasing with probe photon energy above 3.0 
eV as can be inferred from the imaginary part of the GaN 
model dielectric function in Fig. |21 Hence, at shorter 
wavelengths, the probe can not penetrate into the sam- 
ple to detect the phonons. 



E. Coherent control of differential reflectivity 

Coherent optical control of coherent acoustic phonon 
induced differential transmission oscillations by means 
of femtosecond laser sources has been experimentally 
demonstrated in InGaN/GaN multi-quantum well struc- 
tures in Refs.QandjlS. In Ref. 7 control of coherent differ- 
ential transmission oscillations is achieved in a 14 period 



12 



Probe Wavelength (nm) 

450 400 350 300 
2.8 I 




^5 3.0 3.5 4.0 4.5 
Probe Photon Energy (eV) 



FIG. 10: (Color Online) Index of refraction in GaN (upper 
panel) and period of coherent phonon induced reflectivity os- 
cillations (lower panel) as functions of probe photon energy. 
The solid curve in the lower panel is based on Eq. 15511 while 
the dots are obtained from the simulations described in the 
text and presented in Fig. |S| 



InGaN/GaN multi-quantum well sample by means of a 
control pulse whose time delay and intensity could be 
controlled independentlyi^ 

In this paper, we investigate the feasibility of optical 
control of the differential reflectivity oscillations in the 
pin diode structure shown in Fig.^ We consider a simple 
control scenario in which a 3.2 eV pump pulse is followed 
by a delayed 3.2 eV control pulse with the same Gaus- 
sian lineshape as the pump but with possibly a different 
fluence. 

Figs. ^2 (a) and (b) illustrate how the simulated time 
resolved differential reflectivity signals seen in Fig.|Hlcan 
be controlled by means of a delayed control pulse for 
probe photon energies of 3.0 and 3.2 eV, respectively. 
The curves at the bottom of the figures are the time re- 
solved differential reflectivity signals in the absence of 
a control pulse and are the same as the corresponding 
curves in Fig. |S| The upper two curves in Figs. 1111 (a) 
and (b) are simulated transient differential reflectivity 
signals in a pump-probe experiment with a delayed con- 
trol pump in phase and out of phase with the reflectivity 
oscillation period. Both initial and control pulses are 1 
ps Gaussian pulses with linewidths of 4 meV and a pho- 
ton energy of 3.2 eV. We varied the fluence of the control 
pulse relative to the initial pulse in order to maximize 




20 40 60 80 100 20 40 60 80 100 
Time (ps) Time (ps) 



FIG. 11: (Color Online) Differential reflectivity as a function 
of probe delay for probe photon energies of (a) 3.0 eV and (b) 
3.2 eV. The curves are offset for clarity. The pump photon 
energies are 3.2 eV in all cases. The lower curves show the 
differential reflectivity in the absence of a control pulse while 
the middle and upper curves are, respectively, the differential 
reflectivity with a 3.2 eV control pulse in phase and out of 
phase with the reflectivity oscillations at the probe photon 
energy. The best control results for destructive interference 
are achieved with control pulses having approximately twice 
the fluence of the initial pump pulse. 

destructive interference in the differential reflectivity for 
the out of phase cases. We got the best results using 
a control pump having twice the fluence of the initial 
pump. The middle two curves in the Figs. 1111 (a) and 
(b) show the time resolved differential reflectivity when 
the control pulses are in phase with the period of the 
reflectivity oscillations. As we can see the time resolved 
differential reflectivity curves in this case are enhanced 
due to constructive interference. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have simulated time-dependent two- 
color differential reflectivity experiments on a pseudo- 
morphically strained piezoelectric InGaN/GaN multi- 
quantum well in a pin diode including bothe the back- 
ground and oscillating coherent phonon signal. We 
perform simulations in which the pump laser creates 
electron-hole pairs in the strained InGaN wells by pump- 
ing below the GaN band gap. We conclude that very 
large amplitude differential reflectivity oscillations re- 
sult from the generation of coherent acoustic phonon 
wavepackets in the multi-quantum well and their sub- 
sequent propagation into the diode. The propagation 
of these coherent, localized strain pulses into the diode 
results in a position- and frequency-dependent dielec- 
tric function which propagates at the LA phonon sound 
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speed. The diode with the propagating disturbance acts 
Uke a Fabry-Perot interferometer and the period of the 
differential reflectivity oscillations, determined from a 
simple geometrical optics argument, is proportional to 
the probe wavelength. 

To take into account the time dependent background 
differential reflectivity, we modeled the two-color pump- 
probe reflectivity experiments in a Boltzmann equa- 
tion formalism. Electronic structure in the InGaN/GaN 
multi-quantum well was calculated using k ■ p theory. 
We included 1) photogeneration of hot carriers in the 
multi-quantum wells by a pump laser and 2) their sub- 
sequent cooling by emission and absorption of confined 
LO phonons. Nonradiative recombination of electron- 
hole pairs was also included in a simple relaxation time 
approximation. We found that the background differen- 
tial reflectivity signal in our simulated diode for probe 
photon energies near the band gap was smaller than the 



amplitude of the coherent phonon induced differential re- 
flectivity oscillations. 

We demonstrate that coherent optical control of the 
time resolved differential reflectivity can be achieved. We 
have investigated the feasibility of optically controlling 
the differential reflectivity signal using a delayed control 
pulse with the same photon energy and linewidth as the 
pump but with a different fluence. By applying the con- 
trol pulse with a time delay in phase and out of phase 
with the differential reflectivity oscillations, the differ- 
ential reflectivity signal can be selectively amplified or 
suppressed. 
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